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ABSTRACT 

We report the discovery of a new young stellar cluster in the outer Galaxy located at 
the position of an IRAS PSC source that has been previously mis-identified as an exter¬ 
nal galaxy. The cluster is seen in our near-infrared imaging towards IRAS 04186+5143 
and in archive Spitzer images confirming the young stellar nature of the sources de¬ 
tected. There is also evidence of sub-clustering seen in the spatial distributions of 
young stars and of gas and dust. 

Near- and mid-infrared photometry indicates that the stars exhibit colours com¬ 
patible with reddening by interstellar and circumstellar dust and are likely to be low- 
and intermediate-mass YSOs with a large proportion of Class I YSOs. 

Ammonia and CO lines were detected, with the CO emission well centred near 
the position of the richest part of the cluster. The velocity of the CO and NH 3 lines 
indicates that the gas is Galactic and located at a distance of about 5.5 kpc, in the 
outer Galaxy. 

Herschel data of this region characterise the dust environment of this molecular 
cloud core where the young cluster is embedded. We derive masses, luminosities and 
temperatures of the molecular clumps where the young stars reside and discuss their 
evolutionary stages. 

Key words: stars: formation - infrared: stars - sub millimetre: ISM - ISM: clouds - 
ISM: Individual objects: IRAS 04186+5143 - ISM: dust, extinction. 


1 INTRODUCTION 

It has been well established that star formation occurs across 
the Galactic disc and at different Galactocentric distances. 
Both in the inner and in the outer Galaxy, young stellar 
clusters still partly embedded in the dense gas and dust in 
molecular clouds have been found (e.g. Tapia et al. 1991; 
Strom, Strom & Merrill 1993; McCaughrean & Stauffer 
1994; Horner, Lada & Lada 1997; Luhman et al. 1998; San¬ 
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tos et al. 2000). They represent current active star forma¬ 
tion sites. The star formation activity seen throughout the 
Galactic disc is possible due to the relatively large amounts 
of dust that shield young forming stars from the heating of 
the external interstellar radiation Held. The dust produces 
high values of extinction resulting in lines of sight across the 
Galactic disc that are highly opaque in the optical wave¬ 
lengths. 

Among other effects, the large values of dust extinction 
along the Galactic disc make difficult the task of achiev¬ 
ing a complete census of the Milky Way neighbour satellite 
galaxies. Even with our best instruments, we may not have 
found and catalogued correctly all the stellar systems com¬ 
ponents of the Milky Way and its neighbour galaxies. Con¬ 
versely, Galactic stellar systems can be wrongly classified as 
extragalactic neighbours. As an example of this fact, Martin 



2 J. L. Yun, D. Elia, A. A. Djupvik, J. M. Torrelles, and S. Molinari 


et al. (2004) have claimed the discovery of remains of a satel¬ 
lite dwarf galaxy, a claim that was subsquently challenged 
(Momany et al. 2006). More recently, the search for dwarf 
satellite galaxies, both of the Milky Way and of Andromeda 
continues (e.g. Conn et al. 2012; Sesar et al. 2014). 

The detection and characterisation of star formation 
sites in an early stage have strong implications on the struc¬ 
ture and evolution of the Galaxy. However, the study and 
census of star formation sites in the outer Galaxy, and spe¬ 
cially at large distances, has received less attention and cov¬ 
erage when compared to the inner Galaxy and the solar 
neighborhood. 

IRAS 04186+5143 is an IRAS PSC source in the outer 
Galaxy that appears classified in the SIMBAD data base as 
“2MASX J04223304+5150346 - Galaxy”. This means that 
it is listed in 2MASX (the 2Micron All-Sky Survey extended 
source catalogue) as being an extragalactic source. It is also 
an extended subnmi and far-infrared source, having been 
detected in the submm continuum and listed in the SCUBA 
legacy catalogues (Di Francesco et al. 2008), and also de¬ 
tected by the Herschel satellite (Ragan et al. 2012). In addi¬ 
tion, mid-infrared spectral features have been seen towards 
this region (ISOSS J04225+5150 East) using Spitzer (Pitann 
et al. 2011), and Birkmann (2007) has derived a kinematic 
distance of 5.5 kpc. Furthermore, Sunada et al. (2007) found 
no water maser emission in their survey. All the authors 
above clearly refer to this source as a Galactic object. 

As part of our study of young embedded clusters in the 
outer Galaxy (e.g. Yun et al. 2009; Palmeirim & Yun 2010), 
we have conducted observations (near-infrared JHKg imag¬ 
ing, and millimetre CO line) towards IRAS 04186+5143. 
These observations revealed the presence of a young stellar 
population embedded in a molecular cloud core. We report 
here our near-infrared discovery of a young stellar cluster 
seen towards IRAS 04186+5143, and exhibiting evidence of 
sub-clustering. In addition, we use new CO data, as well 
as Herschel observations, archive Spitzer, and archive am¬ 
monia VLA data to characterise the molecular environment 
and the young stellar population. Section 2 describes the 
observations and data reduction. In Section 3, we present 
and discuss the results. A summary is given in Section 4. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 Near-infrared observations 

Near-infrared (J, H , and Ks) images were obtained on 
September 8th 2009 using the Nordic Optical Telescope 
near-IR Camera and Spectrograph (NOTCam). The detec¬ 
tor was the 1024 x 1024 x 18 micron Hawaii science grade 
array (SWIR3). The wide-held camera (0.234"/pix) was 
used, and the observations were performed using a ramp¬ 
sampling readout mode. Every sky position was integrated 
for 36 ( K ) or 48 (J and H ) seconds, reading out the ar¬ 
ray every 6 ( K ) or 8 (J and H) seconds, and using the 
linear regression result of the 6 readouts. The raw images 
were corrected for non-linearity using a pixel-by-pixel cor¬ 
rection model available for NOTCam. Differential twilight 
Hats were used for Hat-fielding. All images were bad-pixel 
corrected, flat-fielded, sky-subtracted, distortion corrected 
(using a model of the WF-camera distortion), shifted, and 


combined to one deep image per filter. The total integra¬ 
tion time in the final images is 648 ( K) and 816 (J and H) 
seconds. 

Point sources were extracted using daofind with a de¬ 
tection threshold of 5a. The images were then examined for 
false detections and a few sources were eliminated by hand. 
Aperture photometry was made with a small aperture (ra¬ 
dius = 3 pix, which is about the measured full-width-at- 
half-maximum of the point spread function) and aperture 
corrections, found from 20 bright and isolated stars in each 
image, were used to correct for the flux lost in the wings of 
the PSF. The error in determining the aperture correction 
was < 0.02 mag in all cases. These errors were added to the 
value MERR, which is created by the IRAF task phot. A 
total of 848, 984, and 786 sources were found to have fluxes 
in J, H, and Ks, respectively, and errors ok s < 0.25 mag. 

We used the 2MASS All-Sky Release Point Source Cat¬ 
alogue (Cutri et al. 2003; Skrutskie et al. 2006) to calibrate 
our observations. The JHKs zeropoints were determined 
using 2MASS stars brighter than Ks = 14.5 mag. The stan¬ 
dard deviations of the offsets between NOTCam and 2MASS 
photometry are 0.04, 0.06, and 0.05 mag, in J, H, and Ks, 
respectively. We estimate the completeness limit of the ob¬ 
servations to be roughly 19.0, 18.5, and 18.0 magnitudes in 
J, H , and Ks, respectively. 

2.2 Millimetre line observations 

The region around the position of the IRAS source was 
mapped using the single-dish OSO 20-m radiotelescope (On- 
sala, Sweden) in 2009 April. Three maps were obtained in 
the rotational lines of 12 CO(1-0), 13 CO(1-0), and CS(2-1) 
at 115.271, 110.201, and 97.981 GHz, respectively. Since the 
telescope half-power beam width (HPBW) is 33”at 115 GHz, 
we decided to obtain the maps with a grid spacing of 30”, 
centered on the IRAS coordinates and composed by 5 x 5 
pointings for the two CO lines, and by 3 x 3 pointings for the 
CS(2-1) line. The typical integration time was 120, 240, and 
300 s for 12 CO(1-0), 13 CO(1-0), and CS(2-1), respectively. 

A high resolution 1600-channel acousto-optical spec¬ 
trometer was used as a back end, with a total bandwith of 40 
MHz and a channel width of 25 KHz that, at the observed 
frequencies, corresponds to 0.065, 0.068, and 0.076 km s _1 , 
respectively. The spectra were generally taken in dual beam 
switching mode, except for the 3x3 innermost portion of 
the 13 CO(1-0) map, observed in frequency switching mode. 
The antenna temperature was calibrated with the standard 
chopper wheel method. Pointing was checked regularly to¬ 
wards known circumstellar SiO masers; pointing accuracy 
was estimated to be 3”rms in azimuth and elevation. 

The data reduction consisted of a typical pipeline for 
mm spectra: first, a folding operation was applied only to 
frequency-switched spectra; then, the baseline has been fit¬ 
ted by a third-order polynomial, and subtracted from the 
spectra (the resulting rms noise per channel is 0.83, 0.35, 
and 0.09 K for 12 CO(1-0), 13 CO(1-0), and CS(2-1), respec¬ 
tively); finally, in all spectra the antenna temperature Ta 
was translated in main beam temperature Tmb dividing by 
the telescope main beam efficiency factor timb- This param¬ 
eter is generally quoted as a constant of the telescope, but 
instead it can vary with the elevation of the source; since 
this variation is evaluated and provided by the OSO 20m 
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telescope system at each pointing we chose to divide 

each spectrum by its peculiar tjmb (i,j) value. 


2.3 VLA ammonia observations 

Simultaneous observations of the NHs(l,l) and NH 3 (2,2) 
lines (rest frequencies 23.694495 GHz and 23.722633 GHz, 
respectively) were carried out with the Very Large Ar¬ 
ray (VLA) of the National Radio Astronomy Observatory 
(NRAO ) 1 in the D configuration during 2003 April 19 
(project AK562; NRAO public archive data). A bandwidth 
of 3.1 MHz with 63 spectral channels of 48.8 kHz width 
(~0.62 km s -1 at A = 1.3 cm) was selected for each am¬ 
monia line. The center channel velocity was set at Vlsr = 
—43.7 km s -1 , covering a total velocity range —62.9 km s -1 
< Vlsr ;$ —24.5 km s -1 . The absolute coordinates of the 
phase center were a(J2000) = 04 h 22 m 34.358 s , <5(J2000) = 
51°50 , 51.0 ,/ , which is ~20” northeast from the nominal po¬ 
sition of IRAS 04186+5143. The observing on-source time 
was ~3 hours. 0542+498 was used as flux calibrator, assum¬ 
ing a flux density of 1.78 Jy at 1.3 cm. The phase calibrators 
were 0359+509 and 4C50.11, with bootstrapped flux densi¬ 
ties 9.04+0.02 Jy and 8.46+0.05 Jy at 1.3 cm, respectively. 
Calibration and imaging was made using the Astronomical 
Image Processing System (AIPS) software of NRAO. The 
resulting synthesised beam size was ~ 3.3” x 2.9” (p.a. = 
78°) with the uv data naturally weighted. An rms per spec¬ 
tral channel of ~ 1.8 mjy beam -1 was obtained in the im¬ 
ages. We estimate that the absolute positions are accurate 
to ~0.5". 


2.4 Spitzer observations 

We searched the Spitzer data archive, and found observa¬ 
tions of the IRAS 04186+5143 region in all the four bands 
(3.6, 4.5, 5.8, and 8 pm) of the InfraRed Array Camera 
(IRAC, Fazio et al. 2004), and at 24 /./in of Multiband Imag¬ 
ing Photometer for Spitzer , being part of the program PID 
20444, executed on 2005 September 20. For IRAC observa¬ 
tions, we used the basic calibrated data (BCD) images pro¬ 
duced by the S18.7.0 pipeline of the Spitzer Science Center: 
96 dithered frames with a 10.4 s exposure, and 32 with 0.4 s 
are available. After having removed the residual muxbleed 
artifacts from the single frames (Hora et al. 2004), we com¬ 
bined into mosaics using the MOPEX software (Makovoz 
& Khan 2005) to obtain two mosaics for each band, corre¬ 
sponding to the long- and to the short-exposure time. The fi¬ 
nal maps have a size of ~ 5'.5x5'.5, and a scale of 0”. 6 /pixel. 

Point source detection and photometry extraction was 
performed with MOPEX as well (Makovoz & Marleau 2005), 
independently at each band. When a source is present both 
in the long- and in the short-exposure image, the photometry 
taken from the latter is considered as more reliable. After 
band merging (based on simple spatial association) a four- 
band catalog of 1020 entries (having at least a detection in 
one of the bands) has been obtained. I 11 particular, sources 
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with detections at the four bands are 215, whereas sources 
detected in bands 2, 3 and 4 are 221. 

2.5 Herschel observations 

IRAS 04186+5143 was observed in the far-IR within the 
Herschel Infrared Galactic Plane Survey (HI-GAL, Moli- 
nari et al. 2010), a Herschel open-time key project which 
mapped the Galactic plane with the Photodetector Array 
Camera and Spectrometer (PACS, 70 and 160 pm; Poglitsch 
et al. 2010) and the Spectral and Photometric Imaging Re¬ 
ceiver (SPIRE, 250, 350 and 500 pm; Griffin et al. 2010) 
instruments on board the Herschel satellite (Pilbratt et al. 
2010). The HI-GAL observations are arranged in tiles of 
~ 2.3° x 2.3° taken at each of the five wavelengths. IRAS 
04186+5143 can be found in the HI-GAL field centred at 
[t,b\ = [152°,+1°] and identified as Field_151_0 in the Her¬ 
schel Science Archive, observed by Herschel on February 
13th 2012 in PACS+SPIRE parallel mode at a scan speed 
of 60 arcsec s -1 . The data were reduced using the UNIMAP 
pipeline, a map maker developed within the Hi-GAL project 
(Piazzo et al. 2015). The maps have pixel sizes 3.2, 4.5, 6 , 
8 and 11.5 arcsec at 70, 160, 250, 350, 500 pm, respectively. 
As in Elia et al. (2013), the astrometry of the maps was 
checked by comparing the positions of several isolated com¬ 
pact sources appearing in both the 70 pm map and in the 
WISE survey (Wright et al. 2010) at 22 pm. Finally, a zero- 
level offset, obtained by comparing the Herschel data with 
Planck and IRAS data, following Bernard et al. (2010), was 
evaluated and added to Herschel maps at each band. 

Compact source extraction and photometry have been 
performed using the Curvature Threshold Extractor pack¬ 
age (CuTEx, Molinari et al. 2011), adopting the same pre¬ 
scriptions and settings used for the general Hi-GAL compact 
source catalog (Molinari et al. 2014, in prep.). The subse¬ 
quent band merging procedure has been carried out based 
on simple spatial association criteria (see, e.g. Elia et al. 
2010 ). 


3 RESULTS AND DISCUSSION 
3.1 The infrared morphology 

Figure 1 presents the NOTCam JHKs near-infrared colour 
composite image obtained towards IRAS 04186+5143. A 
higher concentration of “red” stars (much brighter in the 
A's-band than in the J or i7-bands) is seen close to the 
centre of the image. The ability to resolve most stars in this 
concentration (with possible exceptions at the most crowded 
region) argues in favour of these sources being Galactic. In 
addition, the location of this concentration of stars coin¬ 
cides with the position where the molecular gas, traced by 
CO, peaks (see below), marking the presence of a molecular 
clump. This good spatial coincidence of red stars and molec¬ 
ular gas strongly argues in favour of their association. Thus, 
this image reveals a young stellar cluster still embedded in 
a dense cloud core located in the outer Galaxy. 

Under a closer look, Figure 1 hints at the presence of 
sub-clustering. In fact, the red stars appear to cluster around 
the centre of the frame, but also around a more western 
point. Interestingly, as we show below, the column density 
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Figure 1. J (blue), H (green), and h's (red) colour composite 
image towards IRAS 04186+5143 covering 4' X 4'. North is up 
and East to the left. Note the concentration of red stars towards 
the centre of the image and in the northwest corner. The central 
concentration itself exhibits two groups (sub-clusters) of red stars. 

of the molecular cloud core seems to have a secondary peak 
west of the centre, coincident with the location of the west¬ 
ern red stars. Furthermore, a third smaller group of red stars 
is seen towards the northwest corner if the image. These 
red sources are also seen in the WISE (Wright et al. 2010) 
archive data base colour composite image shown in Figure 2. 

Moving to the far-infrared, Figure 3 shows a Herschel 
RGB (70-160-350) fim composite map of this region at a 
larger scale. The green square indicates the region observed 
by the NOT. IRAS 04186+5143 appears embedded in its en¬ 
vironment, with fainter filaments connecting bright clumps 
and joining at the location where star formation is most ac¬ 
tive (cf. Schneider et al. 2012), a morphology typical of a 
Galactic star forming region. 


3.2 Molecular gas morphology and kinematics 

3.2.1 CO and CS 

In the three maps of the CO transitions, all observed spectra 
clearly show line emission ( S/N > 5). In particular, the 
12 CO(1-0) spectra show a double-peak appearance, that can 
be easily interpreted as self-absorption after comparing with 
the peak positions of 13 CO(1-0) and CS(2-1). In Figure 4, 
the corresponding three spectra observed towards the (0,0) 
position (i.e. the IRAS source location) are overplotted, and 
it looks evident that the 13 CO(1-0) and CS(2-1) peaks lie 
in the range where the 12 CO(1-0) shows a dip between its 
two peaks. Therefore, as it is, this transition cannot be used 
to derive gas physical parameters but, on the other hand, 
indicates a high column density cloud. 

We first used the CS(2-1) line peak to derive the Vlsr 
of the cloud: the center of the Gaussian fit, in the (0,0) 
position, corresponds to Vlsr = —38.2 km s _1 . Similarly, 
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Figure 2. WISE colour composite image centred at 
IRAS 04186+5143 covering 5' x 5'. North is up and East 
to the left. Sub-clusters “E”, “W”, and “NW”, mentioned in the 
text, are marked. “E”, “W”, “C”, and “NW” represent also the 
clumps detected in Herschel’s images (see below). 
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Figure 3. Large-scale Hi-GAL three-colour composite image 
(blue, 70 ^m; green, 160 /xm; red, 350 fim) of the region around 
IRAS 04186+5143. The green box encompasses the area observed 
in the near-infrared. The ruler at the bottom-left corner corre¬ 
sponds to 5 pc (at a distance of 5500 pc). 
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Figure 4. Spectra of the three observed lines toward the centre position of the CO map, corresponding to the location of the IRAS 
source. Blue: 12 CO(1-0); red: 13 CO(1-0); green: CS(2-1). The dashed lines represent the Gaussian fits of the observed lines (see text). 
The vertical grey dotted line indicates the peak position of the Gaussian fit of the CS(2-1) line. 



RA offset (") 

Figure 5. Contour maps (smoothed) of the 12 CO(1-0) and the 13 CO(1-0) integrated intensities (dashed blue and solid red lines, 
respectively), superimposed on the Ag-band image. The OSO 20 m beam at 115 GHz is displayed in the left bottom corner. Labels are 
in units of K km sec -1 . The white crosses indicate the positions observed. 
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for the 13 CO(1-0), we obtained Tlsr = —38.0 km s^ 1 . 
Finally, fitting a Gaussian line profile to the 12 CO(1-0) 
line wings (cfr. Kramer et al. 2004, see Figure 4), yields 
Vlsr = —38.5 km s - . Given the good agreement of all these 
velocities, we adopted the value Vlsr = —38.2 km s -1 . 

According to the circular rotation model by Brand & 
Blitz (1993), this value of Vlsr at the Galactic coordinates 
of this source corresponds to a heliocentric distance dn = 
5.5 kpc, and a Galactocentric distance dc = 13.6 kpc. This 
is in good agreement with the distance quoted by Pitann 
et al. (2011); Birkmann (2007). Thus, the projected sizes of 
the regions mapped in CO (both isotopes) and CS turn out 
to be ~ 4 pc and ~ 2.4 pc, respectively. 

The maps of the integrated intensity / = J Tmb dv of 
12 CO(1-0) and 13 CO(1-0), obtained in the ranges between 
—47 and —27 km s^ 1 and —45 and —30 km s _1 , respectively, 
reflect the different role of these tracers and the presence of 
saturated (self-absorbed) lines in the first one (see Figure 5). 
In fact, the 12 CO(1-0) intensity appears arranged in a single 
“clump” peaked on the (0,0) position, whereas the 13 CO(l- 
0 ) shows a further increase towards the west side of the map, 
corresponding to the second “sub-cluster” (labelled “West”) 
that can be noticed in the infrared images (Figures 1, 2). The 
“NW” sub-cluster lies outside the range of the CO maps. 

3.2.2 Ammonia 

We detected the main component and the inner satel¬ 
lite lines of the NH 3 ( 1 , 1 ) transition in the velocity range 
—39.4 km s -1 < Vlsr iS —36.9 km s -1 , with maximum in¬ 
tensity in the —38.8 km s _1 velocity channel, in very good 
agreement with the CO and CS data. NH3(2,2) emission was 
not detected. In Figure 6 we show the observed NH3(1,1) 
spectrum obtained toward the peak position of the emis¬ 
sion. 

In order to optimise the signal-to-noise ratio of the dis¬ 
tribution of the NH 3 ( 1 , 1 ) emission in the region, we have 
made images with natural weighting and a restoring beam 
of 5”. The corresponding contour maps of different veloc¬ 
ity channels are shown in Figure 7. These contour maps 
reveal that ammonia emission is not detected towards the 
IRAS source position where the main cluster is located. This 
is to be expected as objects embedded in ammonia cores 
are likely to be in earlier stages with little or no near/mid- 
infrared emission detectable. Instead, high-density gas, com¬ 
monly traced by ammonia, is seen here forming structures 
that could represent arcs around the cluster, possibly re¬ 
maining gas from the original parent star-forming core. In 
any case, the detection of ammonia at a Vlsr coincident 
with that of CO corroborates the presence of Galactic dense 
molecular gas. 

3.3 Young stars: the near-infrared view 

Table 1 (full table provided on-line only) gives the photom¬ 
etry of all NOTCam sources detected in the images. 

Figure 8 shows the histogram for the observed (H — Ks ) 
colours of the sources detected in both the H and the Ks- 
band images. The corresponding histogram of a normal star 
field (constituted by main-sequence stars and without the 
presence of embedded clusters) would be approximately a 



Figure 6. NH 3 (1,1) spectrum observed toward the peak position 
of the ammonia emission [a(J2000) = 04 ,l 22 m 33.82 s , 5(J2000) = 
51°51'05.5"]. 
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Figure 7. Velocity channel contour maps of the NH 3 (1,1) main 
line component towards the region of IRAS 04186+5143. The 
LSR velocity is indicated. Contours are -4, 4, 5, 6 , 7, 8 , 9 times 
1.8 mjy beam -1 , the rms noise of the maps. The synthesized 
beam is 5 ". The nominal position of the IRAS source is shown 
by a cross (the size of the cross indicates its position error). 

Gaussian, with the spread in (H — Ks) values about the 
peak value stemming from the range of intrinsic colours of 
main-sequence stars and to low values of variable foreground 
extinction in the lines of sight of each source. However, the 
histogram for the observed (H — Ks) colours of the sources 
towards IRAS 04186+5143 shown in Fig. 8 clearly deviates 
from a Gaussian, exhibiting a red tail (with possibly a second 
peak) representing sources with large values of (H — Ks). 
This observed excess near-IR emission must be due to the 
presence of embedded young stars. 

It is instructive to check the spatial segregation of these 
sources by colour. A Gaussian fit around the peak, exclud¬ 
ing the red tail sources, yields a mean of (H — Ks) = 0.20, 
with a standard deviation of 0.22 for the “blue” peak main- 
sequence field stars. The spatial location of these peak 
“blue” sources is seen in Fig. 9 where they are represented 
by blue open circles. Red filled circles, on the other hand, 
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Table 1. NOT sources towards IRAS 04186+5143 (full table pro¬ 
vided on-line) 


ID 

# 

R.A. 

(°) 

Dec 

(°) 

m j 

m h 

m Ks 

1 

65.58004 

51.83786 

17.45 

16.12 

15.30 

2 

65.58015 

51.84943 

19.41 

18.47 

18.27 

3 

65.58021 

51.85655 

18.66 

17.27 

16.61 

4 

65.58033 

51.83777 

17.68 

16.46 

15.52 

5 

65.58053 

51.85194 

17.62 

16.98 

16.93 

6 

65.58061 

51.85741 

19.58 

18.67 

18.51 

7 

65.58105 

51.82590 

20.26 

19.11 

18.23 

8 

65.58124 

51.86619 

19.79 

18.79 

18.58 

9 

65.58126 

51.86396 



19.10 

10 

65.58178 

51.81469 



19.36 

11 

65.58207 

51.85682 

15.21 

14.82 

14.62 

12 

65.58224 

51.84251 

20.52 

18.62 

17.48 

13 

65.58228 
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Figure 8. Histogram of the observed (H — Kg) colours. The 
well-defined peak is composed of foreground field sources (“blue” 
sources). The red wing of the distribution is composed of sources 
spatially concentrated in the region where the molecular cloud is 
present (see Fig. 9). 


represent sources in the red tail of the ( H — Ks) histogram, 
the “red” sources with ( H — Ks) > 0.75, that is 2.5<r away 
from the blue peak. The blue sources are scattered randomly 
and uniformly across the image, whereas the red sources are 
concentrated in the region of the molecular cloud. Taken to¬ 
gether, these results strongly indicate that most red sources 
are objects associated with the cloud, either located behind 
the cloud (very few sources given the location of the cloud 
in the far outer Galaxy), or being embedded in the cloud 
and possibly containing near-infrared excess emission from 
circumstellar material. 

Blue sources, on the other hand, may be composed 



100 50 0 -50 -100 

A R. A. 


Figure 9. Spatial distribution of all the sources seen both in the 
H and in the Kg band images. Blue open circles represent sources 
with values of (H — K) < 0.75. Red filled circles represent sources 
with values of (H — K) + 0.75. This plot is centred on the IRAS 
point source. 


of a mix of foreground field sources and YSOs in a more 
evolved evolutionary stage. These more evolved young stars, 
if present, could be pre-main-sequence objects or even 
intermediate-mass or massive main-sequence stars formed in 
this cloud, which evolve much faster than their lower-mass 
siblings formed at the same time. Their higher masses would 
also contribute to their being bluer and thus not being told 
apart by red colours. 


Using the point sources detected in all three J, H, and 
/Ts-bands, we plotted the near-infrared colour-colour dia¬ 
gram, (J — H) versus (H — Kg), shown in Fig.10. Most stars 
are located within the reddening band where stars appear 
if they are main-sequence stars reddened according to the 
interstellar extinction law (Rieke & Lebofsky 1985), which 
defines the reddening vector (traced here for Ay = 10). 
Pre-main-sequence YSOs, or massive main-sequence stars 
recently formed in this region, which have had time to clear 
the inner regions of their circumstellar discs, he in this re¬ 
gion as well. Giant stars appear slightly above this band. 
On the other hand, the location of stars to the right of the 
reddening band cannot be the result of interstellar redden¬ 
ing alone. They require the effect of emission by hot dust 
such as that in thick circumstellar discs or envelope molecu¬ 
lar cloud cores. Thus, they are likely to be embedded young 
star objects with infrared excess emission from circumstellar 
material (Adams, Lada & Shu 1987). 

For the sources that lie inside the reddening band, the 
highest value of ( H — Ks ) is about 1.4. Using the mean value 
of (H — Ks) ~ 0.3 for field stars (according to Fig. 8), we 
obtain a colour excess E(H — Ks) = 1.1 due to intra-cloud 
extinction. This value corresponds to a maximum visual ex¬ 
tinction produced by the cloud core, through lines-of-sight 
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Figure 10. Near-infrared colour-colour diagram of the region to¬ 
wards IRAS 04186+5143. The solid line represents the loci of 
unreddened main-sequence stars (Bessell & Brett 1988), while 
long-dashed lines indicate the reddening band. The reddening vec¬ 
tor indicates the direction of the shift produced by extinction by 
dust with standard properties. The location of the vertical dashed 
line, derived from Fig.8, is at (H — Kg) = 0.75. The cross in the 
lower right corner represents a typical error bar. 


where stars can be detected, of Av ~ 17 (Rieke & Lebofsky 
1985). 

The location of the vertical dashed line, derived from 
Fig. 8, is at [H — Ks) = 0.75. The two groups of sources, 
blue sources with ( H — I\s) < 0.75 and red sources with 
( H — Ks) ^ 0.75, are very differently distributed on the 
colour-colour diagram. A large fraction of the red sources 
are located outside and to the right of the reddening band, 
whereas the blue sources mostly occupy the inside of the 
reddening band. Thus, most red sources are likely to be 
YSOs. Given their spatial concentration (Fig. 9), these red 
sources with (H — Ks) *2 0.75 together with a fraction of the 
blue sources in this region, seem to represent a small young 
embedded stellar cluster of about 100 young stars forming 
in the molecular cloud. The actual number of stars in this 
young cluster is likely to be larger for at least three reasons. 
Firstly, we chose a conservative value of (H — Ks) = 0.75 
(2.5a above the mean value of the ( H — Ks) of field main- 
sequence stars). Secondly, there are some young stars in 
a more advanced stage of the star formation process, al¬ 
ready free of circumstellar material and exhibiting blue 
colours, thus not pinpointed by our colour selection crite¬ 
rion. Thirdly, we have considered only stars detected in all 
three JHKs bands. 

We can make a rough estimate of the total mass present 
in the stellar content of this cluster. A first estimate re¬ 
sults from assuming 1 Mg stars, yielding about 100 Mg. 
Not much different values are obtained, e.g. adopting a 
Salpeter Initial Mass Function (IMF) and a reasonable range 
of masses: the result is a total stellar mass of about 140 Mg. 

We derive an upper limit for the masses of the YSOs 



Figure 11. Spitzer 3.6 (blue), 4.5 (green), and 5.8 (red) (im 
colour composite image towards IRAS 04186+5143 covering 
about 5.7’ X 5.7’ with Class I sources (red circles) and Class II 
sources (green squares) overlaid. 

present in this cluster in the following mode. The luminosity 
from the cluster region is dominated by the mid and far- 
infrared flux as measured by IRAS. We estimate this Lfir 
to be about 4.6 x 10‘ ? Lg. Assuming that all this luminosity 
is produced by a single star, this would set an upper limit of 
about 9-10 Mg for any massive star present in this cluster. 
We conclude that the young stellar population present in 
this region is composed of low and intermediate-mass stars. 

3.4 The Spitzer view 

Figure 11 presents the Spitzer bands 1-2-3 colour composite 
image obtained towards IRAS 04186+5143. As expected, a 
clear concentration of stars is seen close to the centre of 
the image. These mid-infrared counterparts of the JHKs 
sources appear quite “red” and support the idea that we are 
dealing with a young stellar cluster located far in the outer 
Galaxy. 

In Figure 12, we present the Spitzer four-band colour- 
colour diagram of this region. Following YSO classification 
criteria, e.g. Allen et al. (2004); Gutermuth et al. (2008, 
2009), we used different coloured-symbols to indicate dif¬ 
ferent types of sources: class I (red circles), class II (green 
squares), and class III and field stars (black dots). Our sam¬ 
ple is restricted to those sources with photometric errors 
less than 0.2 mag in all four IRAC bands. We first select as 
Class I sources those that either have ([4.5] — [5.8] > 1) or 
([4.5] - [5.8] > 0.7 and [3.6] - [4.5] > 0.7), referring here to 
colour indices in magnitudes. From the remaining sources in 
the sample, the Class II sources are those which fulfill the 
three requirements: ([4.5] — [8.0] > 0.5) and ([3.6] — [5.8] > 
0.35) and ([3.6]-[5.8] sj (3.5x (([4.5]-[8.0])-0.5)+0.5). This 
is very close to Gutermuth et al. (2008, 2009) except that 
we have ignored the potential extragalactic contaminants 
which should be of marginal importance in the small area 
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Table 2. Photometric data of Spitzer sources towards IRAS 04186+5143 having NOTCam counterparts (full table provided on-line) 


ID 

# 

R.A. 

(°) 

Dec 

(°) 

m j 

m h 

m Ks 

[3.6] 

(jy) 

Flux 

[4.5] [5.8] 

(Jy) (Jy) 

[8.0] 

(Jy) 

[3.6] - [4.5] 

[5.8] - [8.0] 

YSO 

class 

2 

65.58012 

51.85640 

18.66 

17.27 

16.61 

8.861E-05 

6.169E-05 



0.09 



3 

65.58038 

51.83772 

17.45 

16.12 

15.30 

8.344E-04 

1.130E-03 

1.205E-03 

1.258E-03 

0.81 

0.68 

11 

4 

65.58049 

51.85184 

17.62 

16.98 

16.93 

4.274E-05 

2.621E-05 



-0.05 



5 

65.58195 

51.85678 

15.21 

14.82 

14.62 

3.358E-04 

1.990E-04 

1.330E-04 

4.850E-05 

-0.08 

-0.46 


6 

65.58226 

51.87177 

17.60 

17.03 

16.95 

4.411E-05 

2.499E-05 



-0.13 



7 

65.58245 

51.84241 

20.52 

18.62 

17.48 

1.971E-04 

2.828E-04 

3.685E-04 

4.912E-04 

0.88 

0.95 

1 

8 

65.58312 

51.81478 


18.87 

18.13 

3.614E-05 

2.881E-05 



0.24 



9 

65.58325 

51.81277 

17.63 

16.83 

16.58 

6.500E-05 

4.055E-05 



-0.03 



10 

65.58411 

51.84328 



17.79 

8.834E-04 

1.419E-03 

2.952E-03 

1.952E-03 

1.00 

0.19 

1 

11 

65.58423 

51.84535 

16.92 

16.13 

15.90 

1.154E-04 

7.445E-05 



0.01 



12 

65.58426 

51.87494 

18.98 

17.84 

17.30 

8.294E-05 

8.093E-05 

7.059E-05 

6.755E-05 

0.46 

0.59 

11 

13 

65.58432 

51.87392 

17.74 

16.63 

16.35 

1.045E-04 

8.625E-05 

7.429E-05 

7.952E-05 

0.28 

0.71 

11 

14 

65.58480 

51.87159 

17.35 

16.54 

16.36 

1.871E-04 









[5.8] - [8.0] 

Figure 12. Spitzer colour-colour diagram of the region towards 
IRAS 04186+5143. 


explored here. We find 37 Class I and 48 Class II objects in 
the sample of 215 sources with photometry in all four IRAC 
bands. The large fraction of Class I sources is a clear sign of 
an active and very young region. A typical cluster core has 
a Class II / Class I ratio of 3.7 according to the large survey 
by Gutermuth et al. (2009). The presence of a fair number 
of Class I sources that are usually associated with jets or 
outflows may explain the detection by Pitann et al. (2011) 
of diffuse H 2 , [Si II], [Fe II], and [Ne II] Spitzer spectral lines 
that may indicate the presence of shock-excited gas. 

An inspection of the spatial distribution of the Class I 
and Class II sources (see Fig. 11) reveals that they are all 
predominantly found near the dense clumps detected by 
Herschel (see below), with a clear sub-clustering of at least 
the Class I sources. In Figure 13 we show a simple represen¬ 
tation of their spatial distribution through the histogram of 
all projected separations between sources within the same 
group (Kaas et al. 2004). A homogeneous distribution would 
give a broad gaussian, while clustering shows up as struc¬ 
ture, where peaks relate to the clustering scale. A binsize 
of 8” (0.2 pc) was used to ensure a sufficient population 


within the smallest bin. In order to test the stability and 
the statistical significance of the peaks in the distribution, 
we have varied the binsize in steps from about half to about 
1.5 times this value. The strongest peak gives the approxi¬ 
mate diameter of the most populous group. For the Class I 
sample there is an indication of two separate peaks at small 
source separations. The statistical significance is only about 
1 sigma in this histogram, however, and with larger binsizes 
the two peaks merge to a broad maxima across the range 0.8 
-1.7 pc, which is significant to > 4 sigma. The Class II pop¬ 
ulation has a relatively minor peak at 1.2 pc while it has its 
main peak at 2.2 pc. This clearly shows a stronger clustering 
for Class Is than for Class IIs. The sample of near-IR sources 
with [H — K] > 0.75 magnitudes has two peaks at the same 
location as the Class I sources. Lowering the binsize to 4" 
for this more numerous sample, we fine-tune the locations 
of these peaks to 0.9 pc and 1.7 pc with a high statistical 
significance. In addition, this sample has a broader feature 
at 3 pc reflecting the NW sub-cluster distance to the main 
clusters. This latter population is expected to include both 
the Class I and Class IIs and many more sources not re¬ 
solved/detected with IRAC. The fact that it follows so well 
the small scale structure of the Class Is, suggests that most 
of these sources are likely very young cluster members. 

Table 2 (full table provided on-line only) contains the 
photometry of all Spitzer sources having J, H, or Ks coun¬ 
terparts detected in our NOTCam near-IR images, with in¬ 
dication of YSO classification when available. 


3.5 The Herschel view 

In this section we present and discuss the photometric data 
of Hi-GAL sources detected as described in Sect. 2.4. Fig¬ 
ure 14 presents the Hi-GAL maps of the region surveyed 
in the CO transitions, at 70, 160, 250, 350, and 500 fim, 
respectively. 

Two sources (clumps) were clearly detected at all five 
bands (see panel /) : the eastern one (E) clearly corresponds 
to the CO peak and to the main cluster location, while the 
western one (W) is associated with the second cluster. A 
few further and fainter detections, at only one or two bands, 









10 


J. L. Yun, D. Elia, A. A. Djupvik, J. M. Torrelles, and S. Molinari 


100 
80 
60 

40 E- H 
20 


A 




6 8 



0 2 4 6 8 

Separation (pc) 


Figure 13. The distribution of projected separations between 
sources in the different groups: Class I sources (top), Class II 
sources (middle), NIR sources with [H — I\] > 0.75 (bottom). 
The number of separations per bin is shown versus projected sep¬ 
aration (in pc) assuming a distance of 5.5 kpc, using a bin size of 
8” (~ 0.2 pc). 


represent less reliable sources, with the exception of the one 
located midway between E and W, which we designate by 
source C, and which coincides with a further overdensity 
of sources in the JHKs images. In this case, indeed, the 
source C is clearly visible in the 70 /im and 160 pm maps, 
getting confused with E longward of 160 pm. Given the lim¬ 
ited spectral coverage for C, it is not possible to estimate the 
physical conditions of its envelope through a best-fit proce¬ 
dure based on a modified blackbody model, whereas thanks 
to the availability of SPIRE fluxes this procedure is feasible 
for the E and W clumps (Figure 15). 

Before performing this fit, we note that the beam- 
deconvolved observed sizes 9\ noticeably increase at the 
Herschel 350 and 500 pm bands, so that the fluxes mea¬ 
sured at these wavelengths come from larger volumes of dust 
(Motte et al. 2010; Giannini et al. 2012). Thus, we adopted 
a flux scaling strategy, following that of Elia et al. (2010), 
where we impose F\ = F\ x (#25o/#a), f° r ^ = 350, 500 pm. 
This is based on the assumptions that (i) the source is op¬ 
tically thin at A ^ 250 pm (which is used as reference 
wavelength), (ii) the temperature gradient is weak (Motte & 
Andre 2001), and (iii) the radial density profile is p(r) oc r~ 2 
(i.e. M(r) oc r), 

We then fitted a modified black body to the four fluxes 
from 160 to 500 pm. Since the 70 pm flux generally shows an 
excess due to the proto-stellar content of the clump (Schnei¬ 
der et al. 2012), we keep it as an upper limit to further 
constrain the fit. The modified black body expression is: 

F«,= [l - e~ T AB v {Td)^ , (1) 

where T„ is the observed flux density at the frequency v , 
B v (Td) is the Planck function at the dust temperature Td, 
and 12 is the source solid angle in the sky. The optical depth 
is given by 

tv = {v/vq) 13 , (2) 

where vq = c/Ao is the frequency at which r = 1, and 0 
is the exponent of the power-law dust emissivity at large 


wavelengths. Four free parameters are present in the previ¬ 
ous equations. In order to reduce their number, we imposed 
/3 = 2, as in Elia et al. (2013) (see Sadavoy et al. 2013, for a 
detailed justification of this choice), and 12 to be equal to the 
source area observed at 250 pm. Applying this procedure, 
only T and Ao are left free to change 2 . 

The clump mass is subsequently derived from 

M = (d 2 fl/fc re f )r re f , (3) 

(cf. Pezzuto et al. 2012), where k le f and r re f are the opac¬ 
ity and the optical depth, respectively, estimated at a given 
reference wavelength A re f. Here we chose k re i = 0.1 cm 2 g ~ 4 
at A re f = 250 pm (Hildebrand 1983, which already accounts 
for a gas-to-dust ratio of 100 ), while r re f is obtained from 
Equation 2. 

The physical properties of the clumps E and W are re¬ 
ported in Table 3. We remark that these refer to the volume 
enclosed within the source size observed at 250 pm, which is 
15.5” and 9.7” for the clumps E and W, respectively, well be¬ 
low the CO map grid step. The coordinates given are those 
of the counterparts found at 70 pm (namely the shortest 
wavelength available). 

Since the clump E is well-contained in our emission line 
maps, a comparison with CO-derived masses is possible pro¬ 
vided we consider only the central pointing of the CO maps. 
Two methods can be exploited to compute mass estimates 
from our CO observations, as done, for example, in Yun 
et al. (2009). Here, however, we want to consider the Gaus¬ 
sian profiles fitted to the wings of the self-absorbed observed 
lines as genuine recovered line profiles. 

The first method is based on the empirical linear rela¬ 
tion between column density and the 12 CO(l—0) integrated 
intensity, N(Hz) = Xcolco , with Xco larger in the outer 
Galaxy than in the nearby star formation regions (Xco ~ 
2 x 10 20 cm -2 K -1 km -1 s), given by Nakanishi & Sofue 
(2006): Xco[ cm -2 K 1 km -1 s] = 1.4 x 10 20 exp(7?/ll kpc), 
so that in our case Xco = 4.8 x 10 20 cm ' 2 K -1 km ' 1 s. 
The mass at the central pixel obtained using this method is 
M x = 716 M©. 

The second method assumes local thermal equilibrium 
(LTE) conditions, using 12 CO(l—0) as an optically thick line 
and 13 CO(l—0) as an optically thin line (see, e.g. Pineda, 
Caselli & Goodman 2008). The excitation temperature is 
extracted from the peak main beam temperature of the 
CO( 1—0) (T ex = 31 K). Assuming that excitation tempera¬ 
tures are the same for both lines, the column density of 13 CO 
is calculated through LTE relations, (e.g. equations 6 and 4 
of Brand & Wouterloot (1995), respectively). In order to ob¬ 
tain H 2 column densities, a [H 2 / 13 CO] abundance ratio has 
to be assumed. In the far outer Galaxy it is expected to be 
larger than 5 x 10 5 quoted by Dickman (1978) for local dark 
clouds. Adopting the behaviour of the abundance ratio of 
[ 12 CO/ 13 CO] versus the Galactocentric distance suggested 
by Milam et al. (2005), and assuming a [H 2 / 12 CO] abun¬ 
dance ratio of 1.1 x 10 4 (Frerking, Langer & Wilson 1982), a 
ratio of [H 2 / 13 CO] = 1.1 x 10 6 at R = 13.6 kpc is obtained, 


2 Choosing lower values for 0, e.g. 0 = 1.5, can lead to different 
values of clump masses, as in the case of clumps E and W whose 
masses are found to be smaller by a factor ~ 2 compared to those 
below, but with a much worse y 2 . 
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Figure 14. a, b , c, d , e: Hi-GAL maps of the region around IRAS 04186+5143 at 70, 160, 250, 350, and 500 /zrn, respectively, covering 
the area mapped in CO (magenta dashed box). The ellipses denote the FWHM of the sources extracted by CuTEx at the various 
wavelengths. /: CO(l-O) integrated intensity in grey scale, with the sources identified in panels a-b superposed (using the same color- 
wavelength convention of previous panels). Solid ellipses identify the counterparts of the three sources having a reliable SED (see text), 
indicated with “E”, “C”, and “W”, respectively, whereas dashed ellipses identify sources detected in only one or two bands, and discarded 
for further analysis. 


E W NW 



l(«m) m) A(wm) 

Figure 15. Mid and far-infrared spectral energy distributions of the three main Hi-GAL sources found in the investigated region, the 
source “E” being coincident with the location of IRAS 04186+5143. WISE fluxes are represented with triangles, while Herschel ones are 
marked with circles. Fluxes taken into account for the modified grey body fit (shown with a solid curve) are plotted with filled symbols. 
The grey-shaded area represents the integral computed to determine the bolometric luminosity Lboi- 


Table 3. Photometric and physical properties of the main Hi-GAL sources 


Designation 

a 

deg 

6 

deg 

F70 

Jy 

^160 

Jy 

-^250 

Jy 

^350 

Jy 

-F500 

Jy 

0250 

// 

M 

M 0 

T 

K 

Ao 

^m 

-^bol 

L© 

E 

65.6354 

51.8417 

54.4 

51.2 

45.9 

21.4 

9.0 

15.5 

719 

17.1 

83.9 

3150 

C 

65.6252 

51.8435 

10.0 

18.1 








> 510 

w 

65.6105 

51.8454 

6.6 

17.4 

21.9 

10.3 

4.2 

9.7 

435 

15.7 

104.0 

640 

NW 

65.6000 

51.8758 

0.3 

4.6 

8.2 

5.2 

2.3 

25.8 

416 

12.4 

38 

125 
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resulting in a value of Mlte = 635 Mg for the mass derived 
with this method. 

Our Hi-GAL-based mass estimate for the clump E is 
Mh = 719 Mg (see Table 3), similar to the aforementioned 
CO-derived masses. Also, as these three mass estimates have 
the same distance dependence, the agreement among them 
in not affected by the uncertainty in this parameter. In¬ 
stead, the choice of the Xco and [H 2 / 13 CO] implies crucial 
assumptions: the comparison between the CO- and Hi-GAL- 
derived masses should be useful to better calibrate these pa¬ 
rameters, which in this case, would require relatively small 
adjustments. However the k le t value adopted in Equation 3 
is typically used for the inner Galaxy, but the gas-to-dust ra¬ 
tio is expected to increase with decreasing metallicity (con¬ 
ditions that are to be found at large Galactocentric radii), 
leading to large uncertainties on the final gas masses (see 
Mookerjea et al. 2007, and references therein). Thus, adopt¬ 
ing a larger gas-to-dust ratio in this case would imply, ac¬ 
cordingly, a rescaling of Xco and [H 2 / 13 CO]. 

Once a good agreement between mass estimates ob¬ 
tained from the continuum and line maps has been ascer¬ 
tained, we can discuss the relationship between mass and 
size, using the quantities in Table 3. At a distance of 5.5 kpc, 
the angular extent of the clump E corresponds to a phys¬ 
ical diameter of ~ 0.4 pc. This implies a surface density 
E = 1.1 g end 2 , a value that exceeds the theoretical thresh¬ 
old calculated by Krumholz & McKee (2008) for a star form¬ 
ing cloud to be able to form massive (M > 10 Mq) stars. 
This is an interesting finding, testifying the presence of con¬ 
ditions for high mass star formation also in the outer Galaxy, 
and in relatively isolated regions far from giant star forming 
clouds. For the clump W, the surface density is found to be 
even larger, E = 1.7 g cm -2 . 

The relation between the mass of the envelope (i. e. the 
mass we derive from Herschel’s data) and the bolometric lu¬ 
minosity can be used to diagnose the evolutionary stage of 
a far-infrared source (e.g., Molinari et al. 2008; Ma, Tan & 
Barnes 2013; Elia et al. 2013). The L\ m \/M ratio, in partic¬ 
ular, is a distance-independent quantity which is expected 
to rapidly increase during the accretion phase of star forma¬ 
tion. 

Here the bolometric luminosity has been obtained as the 
trapezium-like integral of the observed SED at A < 160 /tm 
(including the fluxes at the WISE bands at 12 and 22 /j,m), 
added to the integral of the best-fitting modified black body 
at A ^ 160 a m, namely the grey-shaded areas represented 
in Figure 15. The values obtained for the clumps E and 
W are reported in Table 3. The Lboi/M ratio amounts to 
4.4 Lq/Mq for clump E, and to ~ 1.5 Lq/Mq for clump 
W. These values, corresponding to the formation of a young 
cluster, cannot be directly compared with models elaborated 
to describe the formation of a single young stellar object 
(Molinari et al. 2008). However, the direct comparison be¬ 
tween clumps E and W shows that the former is likely to be 
in a more evolved stage. Also, the ratio between the bolomet¬ 
ric luminosity and its submillimeter portion Lsubmm derived 
from fluxes longward of 350 /mi, has been used as a further 
evolutionary indicator, expected to be larger at more evolved 
star formation stages (e.g. Andre, Ward-Thompson & Bar- 
sony 1993). The values of this ratio, for the E and the W 
clumps, is found to be 67 and 28, respectively, again confirm¬ 
ing the previous indication. This view is further supported 


by the fact that: (j) clump E has a smaller [70-160] color 
index than clump W; (ii) clump E appears to be slightly 
warmer than clump W (and the temperature estimate does 
not depend on the 70 /mi flux); and (Hi) clump E is less 
dense than clump W and becoming optically thin at shorter 
wavelengths (which could indicate that a larger fraction of 
the gas and dust envelope has already been transferred onto 
the forming stars, or dissipated by proto-stellar activity). 

In Figure 3 it can be noticed that the NW clump (cf. 
Figure 2) is also partially covered by our JHKs maps (the 
red sources at the top right corner of Figure 1 are spa¬ 
tially associated with it) and fully observed by Herschel 
(clearly detected at all bands, with no duplicities), so that 
we can derive the physical properties of its dust and gas 
envelope. A comparison with CO observations, however, is 
not possible since this source lies outside the area mapped 
at OSO. The SED of NW is shown in Figure 15 and the 
results of the fit procedure are reported in Table 3, respec¬ 
tively. This clump appears remarkably fainter and less dense 
(E = 0.2 g cm -2 ) than the others. At the same time, it could 
be going through an earlier evolutionary stage, as testified 
by its low temperature (T = 12.4 K) and luminosity/mass 
ratio (Lboi/M = 0.3 Lq/Mq), and by L ho \/L snhmm = 10. 


4 SUMMARY AND CONCLUSIONS 

• Infrared (JHKs and Spitzer) images of the region to¬ 
wards IRAS 04186+5143 reveal a concentration of stars 
compatible with the presence of a young stellar cluster. 

• The cluster is embedded in a molecular cloud core de¬ 
tected through CO, CS, and NH$ line emission. Our CO 
map reveals the existence of sub-clustering corroborated by 
the spatial distribution of the young stars. 

• At 5.5 kpc (heliocentric distance) in the outer Galaxy, 
and at a galactocentric distance of 13.6 kpc, this population 
of YSOs may be composed by low and intermediate-mass 
stars with a large fraction of Class I sources, a clear sign of 
a young star formation region. 

• Herschel data clearly identify dust clumps coinciding 
with the positions of the sub-clusters. The Herschel-derived 
masses of the main clumps are 719 and 435 Mq, consistent 
with CO -derived estimates. 

• The L/M ratio of the clumps could indicate that the 
larger (E) clump, hosting a larger fraction of the YSOs seen 
in the near-infrared images, is in a more evolved stage of the 
star formation process, having converted more gas into stars 
than the smaller (W) clump. 

• IRAS 04186+5143 is a young stellar cluster forming in 
the outer Galaxy, and not an external galaxy as identified 
in the 2MASS extended source catalog and indicated in the 
SIMBAD data base. 

• A table is provided giving the photometry of all NOT- 
Cam sources detected at least in one near-IR band. 

• An additional table is provided cross-correlating Spitzer 
and NOTCam sources. It contains the mid-infrared photom¬ 
etry of all Spitzer sources with NOTCam J, H, or Ks coun¬ 
terparts. 
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